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ABSTRACT 



> . 

^sO ! A candidate extrasolar planet companion to the young brown dwarf 2MASSW 

^ ! Jl 207334-393254 (2M1207) was recently discovered by Chauvin et al. They find 

that 2M1207 B's temperature and luminosity are consistent with being a young, 
\ ~5 Mj up planet. The 2M1207 system is purported to be a member of the TW 

Hya association (TWA), and situated ~70pc away. Using a revised space mo- 
tion vector for TWA, and improved proper motion for 2M1207, I use the moving 
cluster method to estimate the distance to the 2M1207 system and other TWA 
members. The derived distance for 2M1207 (53 ± 6 pc) forces the brown dwarf 
and planet to be half as luminous as previously thought. The inferred masses 
for 2M 1207 A and B decrease to ~21 M Jup and ~3-4 M Jup , respectively, with 
the mass of B being well below the observed tip of the planetary mass function 
and the theoretical deuterium-burning limit. After removing probable Lower 
Centaurus-Crux (LCC) members from the TWA sample, as well as the prob- 
able non-member TWA 22, the remaining TWA members are found to have 
distances of 49 ± 3 (s.e.m.) ± 12 (la) pc, and an internal ID velocity dispersion of 
O.8I02 kms" 1 . There is weak evidence that the TWA is expanding, and the data 
are consistent with a lower limit on the expansion age of 10 Myr (95% confidence). 

Subject headings: open clusters and associations: individual (TW Hya associ- 
ation) — stars: brown dwarfs — stars: distances — stars: kinematics — stars: 
planetary systems — stars: individual (2MASSW J1207334-393254) 
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1. INTRODUCTION 

Indirect detection methods for finding extrasolar planets have yielded in excess of 150 
candidate planets over the past decade (Marcy et al. 2005). None have been directly im- 
aged, however two have had their light detected through transit studies with the Spitzer 
Space Telescope (HD 209458b and TRes-lb; Charbonneau et al. 2005; Deming et al. 2005). 
Recently, an object has been imaged, and resolved, whose properties appear somewhat con- 
sistent with being a young extrasolar planet: the companion to 2MASSW J1207334-393254 
(2M1207; Chauvin et al. 2004). These objects appear to represent the opening chapter in 
humanity's quest to study the atmospheres of planets beyond our solar system. 

As one of the nearest and youngest brown dwarfs yet identified, 2M1207 has received 
considerable attention since its discovery. Gizis (2002) discovered 2M1207 in a spectroscopic 
survey of red 2MASS sources, and claimed that the object was a ~10 Myr-old, ~25Mj np 
member of the nearby (~55 pc) TW Hya association (TWA; Webb et al. 1999b). Further 
observations of its radial velocity (Mohanty et al. 2003) and proper motion (Scholz et al. 
2005) are roughly consistent with TWA membership. With low resolution spectroscopy, 
Gizis (2002) found 2M1207 to show signs of low surface gravity and strong Ha emission 
(EW = 300A). In their echelle spectroscopy survey, Mohanty et al. (2003) found the Ha 
emission line to be broad and asymmetric, and accompanied by several other Balmer and 
He I emission lines. Mohanty et al. (2003) hypothesize that the brown dwarf is probably still 
accreting from a circumstellar disk. 

While astrophysically interesting in its own right as a representative of the new class of 
young, accreting brown dwarfs (e.g. Muzerolle et al. 2005; Mohanty et al. 2005), it appears 
that 2M1207 may become most famous for being the host "Sun" for the first imaged extra- 
solar planet - if indeed 2M1207B can be called a "planet". Chauvin et al. (2004) discovered 
a faint companion to 2M1207, which has near-IR photometry and a low signal-to- noise-ratio 
spectrum consistent with having a late-L spectral type. Recently, Chauvin et al. (2005) 
and Schneider et al. (in prep.) confirmed that the companion B is indeed co-moving with 
2M1207 A. Debate on the origin and classification of this object is in its infancy. To help 
better constrain the physical nature of this object, I present an improved distance estimate 
to the 2M1207 system through the moving cluster method. The new distance provides more 
accurate luminosities (and inferred masses) for the components of this interesting substellar 
binary. 
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2. ANALYSIS 

Although a trigonometric parallax is not yet available for 2M1207, one can exploit the 
star's putative membership in the TW Hya association to derive a distance using the cluster 
parallax (or "moving-cluster") method (e.g. Atanasijevic 1971; de Bruijne 1999a). With an 
observed proper motion and radial velocity (as well as other supporting evidence), I test 
whether the star is consistent with being a TWA group member. To exploit this technique, 
one needs to take the following steps: (1) estimate the space motion vector for the TWA, (2) 
test whether the observations for 2M1207 (proper motion, radial velocity) are consistent with 
the TWA motion vector, and (3) use the moving cluster method to estimate the parallax from 
the proper motion and TWA space motion data. I address these steps in order. Although 
the rest of the TWA membership is not the focus of this study, I will briefly mention relevant 
results for these systems throughout this analysis. I will also examine whether the expansion 
of the TWA is detectable, and whether it can help constrain the age of 2M1207 and the rest 
of the association. 



2.1. Sample 

The initial pool of candidate TWA members considered in this study are listed in Table 
1. I add to TWA numbers 1 through 25 (Zuckerman & Song 2004) the three, new, low-mass 
candidate members 2M1207 and 2MASSW Jl 13951 1-315921 from Gizis (2002), and SSSPM 
J1102-3431 from Scholz et al. (2005). The TWA members from Webb et al. (1999b, TWA 
1-11) and Sterzik et al. (1999, TWA 12, 13) comprise what I will tentatively call the "classic" 
membership of the TW Hya association. These are young stars which were mostly selected 
due to infrared or X-ray excesses within the immediate vicinity of TW Hya. There has been 
debate regarding the membership for TWA 14-19 (Mamajek & Feigelson 2001; Lawson & 
Crause 2005), since their positions, proper motions, and rotational properties are at variance 
with TWA 1-13. Mamajek & Feigelson (2001) and Lawson & Crause (2005) have suggested 
that TWA 14-19 are probably members of the more distant (~120pc; de Zeeuw et al. 1999) 
and older (~16 Myr; Mamajek et al. 2002) Lower Centaurus Crux OB association (LCC). 
TWA 20 was claimed to be a TWA member by Reid (2003, R03), but rejected by Zuckerman 
& Song (2004) due to its weak Li. As the Li data is not published, and the similarity in 
proper motion between TWA 20 and the other TWA members is quite striking, I retain 
TWA 20 in the candidate pool. TWA 21 through 25 were selected by Song et al. (2003, 
SZB03) due to their strong Li and Ha emission. From Fig. 6 of SZB03, it appears that 
TWA 23 and 25 have positions and proper motions very close to those of TWA 1-13, but 
TWA 21 and 22 are spatially isolated, and TWA 24 has a small proper motion, similar to 
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LCC members. Hence it is not obvious that many of the TWA 14-25 stars were born in the 
same star-formation event as TWA 1-13. I conservatively include only the classic members 
(TWA 1-13) in the initial calculations for estimating the convergent point and space motion 
vector for the TW Hya association. 

2.2. Astrometric Data 

The adopted proper motion and radial velocity data for proposed TWA members, and 
their associated references, are presented in Table 1. I searched the literature and on-line 
catalogs 2 to find the best values for the proper motions and radial velocities of TWA mem- 
bers. To mitigate against the effects of short-term astrometric perturbations by short-period 
companions, I preferentially adopted the long-baseline proper motion with the smallest error 
bars (usually Tycho-2 or UCAC2; H0g et al. 2000; Zacharias et al. 2004) over Hipparcos 
values (Perryman et al. 1997), when available. In a few instances, I calculated new proper 
motions using published positions. I calculated weighted mean radial velocities when mul- 
tiple values were available, or adopted systemic velocities for spectroscopic binaries, when 
available. 

2M1207 has two published proper motion estimates in the literature (Gizis 2002; Scholz 
et al. 2005). The Gizis (2002) proper motion (/i a *,/i«5 = -100,-30 mas yr" 1 ) does not have 
error bars and is based only on a few plate images in the USNO image archive. Scholz et 
al. (2005) estimated a proper motion for 2M1207 of =-78 ± 11 masyr -1 , fi$ = -24 ±9 
masyr -1 . This proper motion estimate included a Chandra pointing, rather than an actual 
measured position, and so is invalid. Omitting the Chandra pointing, R.-D. Scholz has 
calculated a revised proper motion of =-67±7masyr" 1 , /i,5 = -28±ll masyr -1 using 
a least-squares fit with equal weighting (R.-D. Scholz, personal communication). 

As there are large differences in the accuracy between the SuperCOSMOS, 2MASS, and 
DENIS positions (~60mas vs. ~500mas), I recalculated the proper motion using weighting 
by the inverse of the square of the positional errors, following the method of Corbin (1977) 3 . 
The SuperCOSMOS and 2MASS positions are tied to the International Celestial Reference 
System (ICRS) via the Tycho-2 catalog, and so for our purposes they are on the same system. 
In order to estimate a positional error for the SuperCOSMOS positions, I performed a least- 



2 ADS (http://adsabs.harvard.edu/), Vizier (http://vizier.u-strasbg.fr/viz-bin/VizieR), and SIMBAD 
(http: / /simbad. harvard.edu/sim- fid. pi) 

3 The formulae are given in the on-line documentation for the AC2000.2 catalog (Urban et al. 1998) at 
http://ad.usno.navy.mil/ac/ 
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squares fit to the 4 SuperCOSMOS points, and found their scatter consistent with positional 
errors of a m = 143 mas and as = 196 mas. These errors are very consistent with the 
SuperCOSMOS positional errors quoted by Hambly et al. (2001). I corrected the DENIS 
position for the 2MASS-DENIS offset found by Cabrera-Lavers & Garzon (2003), since the 
2MASS positional errors are much smaller than DENIS's (2MASS is tied to the ICRS via the 
Tycho-2 catalog to an accuracy of ~80 mas; Cutri et al. 2003). For the DENIS positional 
errors, I adopted the square root of the 2MASS-DENIS rms differences added in quadrature 
with the 2MASS-ICRS rms residuals (~80mas), giving a aif = 430 mas, and as = 320 mas. I 
estimate the proper motion of 2M1207 to be n aJf =-72 ± 7masyr _1 , /i<5 = -22±9 masyr -1 , 
which is within la of both of Scholz's estimates. The change in the position of 2M1207 over 
time is plotted in Fig. 1. 

2.3. The Space Motion of the TWA 

In order to calculate a cluster parallax for 2M1207, we require an accurate convergent 
point solution for the TW Hya association, to which 2M1207 is proposed to be a member. 
Mean space motion vectors and/or convergent point solutions for the TWA were previously 
estimated by Frink (2001), Makarov & Fabricius (2001, MF01), R03, and SZB03. Considering 
the increase in proposed association membership (SZB03), wealth of new proper motion data 
(UCAC2; Zacharias et al. 2004) and radial velocity data (Torres et al. 2003) made available 
since R03, I will briefly discuss and reanalyze the kinematics of TWA. To estimate the space 
motion for the TWA, I will combine information from 2 different methods: using what little 
data there is regarding the 3D space motion vectors for individual members, as well as 
applying the convergent point method on the classical membership. 

Only four classic TWA members have sufficient data to reliably calculate the 3D space 
motion vector (TWA 1, 4, 9, and 11), and these individual determinations have modest 
errors in any given velocity component (a ~ l-2kms _1 ; Mamajek et al. 2000). Three of 
the systems are binaries, but their systemic velocities are probably accurate to ~1 kms -1 , or 
better. The mean barycentric Galactic space motion vector for these four systems (U,V,W 
= -10.2, -17.1, -5.1 kms -1 ) provides the best estimate of the centroid velocity vector for the 
TW Hya association. 

For helping refine the vertex estimate for the TWA, I will use the convergent point 
method on the classical membership. The convergent point, as calculated only from the 
proper motion data, will also become important when the question of association expansion is 
addressed (§2.6). I approximately follow the convergent point grid technique of Jones (1971). 
In this implementation, I alter Jones's definition of t 2 , following de Bruijne (1999a), and 
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include an intrinsic velocity dispersion term (a v = 1 kms" 1 ), and assumed distance (50 pc) 
in the definition for t 2 (the method is rather insensitive to both input values). Over the 
entire hemisphere a G 0° — 180°, I calculate the t 2 statistic at every 0°.l grid step, and find 
the celestial position which gives the minimum t 2 value. For every grid point, the method 
assumes that this position is the convergent point for the group, and rotates the stellar proper 
motion components (in and n$) to the proper motion directed toward the convergent 
point (fjL v ) and perpendicular to the great circle joining the star and test convergent point 
(fi T ). The method iteratively searches for which test convergent point minimizes the r 
components of proper motion for the input sample. Jones's and de Bruijne's t 2 value can be 
treated statistically as the classic x 2 (Bevington & Robinson 1992). In its iterative search 
for the group convergent point, the method will reject stars contributing the most to the t 2 
statistic, until the position with lowest t 2 value corresponds to a sufficiently high enough x 2 
probability that the best convergent point can not be statistically rejected. For a statistical 
rejection threshold, I adopt a 5% level of significance (i.e. 5% probability of falsely rejecting 
the null hypothesis) following Trumpler & Weaver (1953). 

The TWA stars are sufficiently convergent, and the proper motion errors for the faint 
members are large enough, that a convergent point can be determined for all the classic 
TWA members (#1-13) with a low xl (x 2 l v — 15.9/13; x 2 probability = 25%). For internal 
velocity dispersions of a v > 0.6 kms -1 , the method is able to find a convergent point with 
X 2 probability of >5% without rejecting any of the classic members. If a v = 0.6 kms -1 is 
adopted, TWA 6 (which contributes the most to the t 2 statistic) is rejected, and a sound 
solution is found with the other nuclear members (x 2 probability = 21%). The internal 
velocity dispersion is probably near a v ~ 1 kms -1 , and with this adopted velocity dispersion, 
33% of the classical members contribute At 2 > 1 (similar to how MF01 estimate the velocity 
dispersion). Hence, there is no good reason to remove TWA 6 in the hunt for a statistically 
satisfactory convergent point for TWA 1-13. I will determine a more refined estimate of the 
velocity dispersion for TWA in §2.5. The ability of the technique to give a statistically sound 
convergent point solution (xl — 1) with a v = l kms" 1 already suggests that the velocity 
dispersion of TWA is similar to that of nearby OB associations (Madsen et al. 2002). 

In Fig. 2, I plot the convergent points for subsamples of the TW Hya association, 
as well as previous determinations from the literature. I also plot the convergent points 
for subsamples of the TWA 14-25 membership in Fig. 2. The confidence regions of these 
subsamples are roughly twice as large as than that for TWA 1-13, but not wholly inconsistent 
given the large error bars. Much of the positional deviance of these subsample convergent 
points is owed to TWA 22, which may not be a kinematic TWA member (§2.5). From Fig. 
2, one can conclude that the convergent point for TWA 1-13 (within the dashed confidence 
regions) agrees well with that predicted by the TWA space motion vectors found by R03 and 
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SZB03. The TWA vertex found by MF01 is just outside of the 95% confidence region, and 
seems to be deviant when compared to the values from R03, SZB03, F01, and the results of 
this convergent point analysis. This fact, combined with the finding that most of the stars 
in the MF01 convergent point analysis are not pre-MS (Song et al. 2002), suggests that their 
convergent point and dynamical age for the TWA (8.3 Myr) are not valid. I will discuss the 
expansion age further in §2.6. 

After considering the agreement between the TWA 1-13 convergent point and that 
inferred from the mean TWA space motions from R03 and SZB03, I adopted the following 
fiducial TWA parameters. For the group vertex, I took the weighted mean of the vertices 
from my convergent point analysis of TWA 1-13 (a, 5 = 100° ±10°, -28° ±4°), and the 
individual vertices inferred from the space motion vectors for TWA 1, 4, 9, and 11 (using 
eqn. 10 of de Bruijne 1999a). This analysis assumes zero association expansion, and assumes 
that there is no significant offset between spectroscopic and physical radial velocities - both 
of which are acceptable assumptions at this level of accuracy. The best estimate of the 
convergent point for the TWA is calculated to be (a = 103°.2 ± 1°.5, 5 = -30°.7±l°.5). For 
the mean speed of the classic TWA membership, I adopt the weighted mean barycentric 
speed (v = 21.3 ± 1.3 (s.e.m.) kms -1 ) of TWA 1, 4, 9, and 11, using their astrometry and 
radial velocities in Table 1, and weighted mean Hipparcos and Tycho parallaxes. 



2.4. Is 2M1207 a TWA Member? 

Given the adopted convergent point solution for the "classic" TWA members, and a 
proper motion for 2M1207, one can estimate a membership probability and predict the 
star's radial velocity. Using the updated proper motion for 2M1207 (§2.2), I find that most 
of the motion is indeed pointed towards the convergent point (jj, v = 75 ± 7 mas yr _1 ) and very 
little of it is in the perpendicular direction (/i r = 2 ± 8masyr _1 ). Using the membership 
probability equation from de Bruijne (1999a, his eqn. 23), and adopting a mean cluster 
distance of 50 pc and velocity dispersion of lkms -1 , I estimate a membership probability 
of 98%. This membership probability should be interpreted as: given the proper motion 
errors, 98% of bona fide TWA members are expected to have /i r values more deviant than 
2M1207. That is, the proper motion of 2M1207 is consistent with the null hypothesis (/i r = 0) 
for an "ideal" member. One can also use the predicted and observed radial velocity as a 
check of the moving cluster method. Assuming parallel motion among group members, the 
method predicts the radial velocity as v ra d = v cosA, where v is the speed of the group, and 
A is the angular separation (62°.9±1°.5) between 2M1207 and the convergent point. The 
predicted radial velocity for 2M1207 (+9.7 ± 1.6 km s" 1 ) is within 0.6a of the observed radial 
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velocity measured by Mohanty et al. (2003, +11.2 ± 2kms~ 1 ). Both the proper motion and 
radial velocity data for 2M1207 are quantitatively consistent with TWA membership, and 
its evidence of membership is as strong as that for most of the classical members. 

2.5. Distances 

2.5.1. The Distance to 2M 1207 

If a star belongs to a moving group, its proper motion can be used to estimate its 
distance. The star's moving cluster parallax {w) is calculated as w = A /i v /v sin A, where [i V) 
v, and A are as described before, and A (= 4.74047) is the astronomical unit expressed in 
the convenient units of kmyrs^ 1 (de Bruijne 1999a). Using the values (and uncertainties) 
for fi v , v, and A as given in §2.2 and §2.3, I calculate a cluster parallax for 2M1207 of 
zu = 18.8 ± 2.3 mas, or a corresponding distance of d = 53 ± 6 pc. The only published distance 
estimates to 2M1207 are ~70pc (Chauvin et al. 2004) and 70±20pc (Chauvin et al. 2005). 
Both are photometric distance estimates which force 2M1207A to be an unreddened M8 
star on a 10 Myr-old isochrone. Considering the variations between published evolutionary 
tracks, especially for stars which are young and low-mass (Hillenbrand & White 2004), the 
cluster parallax distance should be considered an improved estimate. 

2.5.2. Distances to TWA Objects: Implications and Final Membership 

The agreement between the trigonometric parallaxes for TWA 1, 4, 9, and 11, and their 
cluster parallaxes are excellent, as shown in Table 2. All the parallaxes are within la of each 
other, with an insignificant weighted-mean zero-point offset of -0.8 ± 1.2 mas, in the sense 
"cluster minus trigonometric". Cluster parallax distances for all TWA member candidates 
are given in column 9 of Table 1 . 

There is a small caveat regarding the cluster parallax distances in Table 1 and Fig. 3 that 
is worth elaborating upon. There have been suggestions that some TWA stars may actually 
be background members of the ~16-Myr-old Lower Centaurus-Crux (LCC) OB subgroup at 
d ~ 110 pc (e.g. Mamajek & Feigelson 2001; Mamajek et al. 2002). The space motion vectors 
for TWA and LCC are very similar, and within roughly ~5kms _1 (Mamajek et al. 2000). 
If one calculates cluster parallax distances to "TWA" objects using the space motion vector 
of LCC (Madsen et al. 2002), the mean distances in column 7 of Table 1 change by less 
than ±5% (rms). This is smaller than the quoted distance errors (typically ~11%). Hence, 
any conclusions based on the distribution of cluster parallax distances (i.e. Fig. 3) are very 
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insensitive to whether individual "TWA" objects are co-moving with either TWA or LCC. 

I plot the cluster parallax distances versus Right Ascension in Fig. 3. Fig. 3 illustrates 
that there appears to be a gap in the distances between LCC members and "classic" TWA 
members near d = 85 pc, effectively splitting the groups spatially. Hence, TWA 12, 17, 18, 19, 
and 24 have distances more consistent with LCC than the other TWA members. Previous 
investigators (MF01, R03) have suggested that the TWA members are clustered at distances 
of ~70 pc, however Fig. 3 suggests that what is really being seen is two detached populations 
of young stars: one at ~50pc (TWA) and one at ~110pc (LCC). The agreement between 
the observed and predicted radial velocities for TWA 12, 17, 18, 19, and 24, are probably 
due to the similarity in space motion between LCC and TWA (see Fig. 2). As it is often not 
clear how "TWA" candidate members have been retained (or rejected) in past studies, there 
may be an observational bias present for the radial velocities of these more distant objects 
to agree well with that of the foreground members. 

As Fig. 3 suggests that some of the "TWA" stars may be more distant members of LCC, 
it is worth reexamining the vertex of the remaining TWA members, including the new brown 
dwarf members TWA 26-28. If the convergent point method (§2.3) is run on the remaining 
members (again assuming a mean distance of 50 pc and o^, = 1 kms -1 ), a somewhat poor 
vertex solution is found (x 2 /^ — 39.8/23; x 2 probability = 1.6%). The biggest contributer 
to the x 2 (contributing a third of the quantity) is the closest TWA candidate - TWA 22. 
If TWA 22 is dropped, the convergent point method shifts by a few a in position, and a 
much more statistically sound vertex is found: a = 100°. 5 ± 5°.0, 5 = —27°. 9 ± 2°. 3, 
X 2 /v = 17.5/22; x 2 probability = 74%. Rejecting further members has negligible effect on 
the vertex, and only pushes the x 2 probability to absurdly higher levels. This remarkable 
reduction in % 2 , upon removal of TWA 22 from the sample, suggests that TWA 22 should 
probably be excluded as a TWA member. Clearly it is a nearby, young star, however it 
does not appear to be a kinematic TWA member. In the initial calculation of membership 
probabilities (column 8 of Table 1; §2.4), TWA 22 had P = 2%- by far the lowest. This 
new a posteriori convergent point estimate is currently the best that can be done purely 
geometrically, i.e. with proper motions alone. It is in excellent agreement with the original 
TWA 1-13 vertex determination, and with the individual vertices for TWA 1, 4, 9, and 11. 

With the sample of "final" TWA members (denoted "Y" or "Y?" in column 11 of Table 
1), one can independently estimate the velocity dispersion a v of TWA based on how well the 
proper motions determine the convergent point. Considering the range of x 2 values for an 
acceptable fit (see discussion in Gould 2003), the final estimate of the velocity dispersion of 
TWA, from the proper motion data alone, is a v = 0.8jl ;| kms -1 . By adopting a v = 0.8 km s -1 , 
the uncertainties on the proper motion-determined convergent point decrease to a a = 4°.2 



-10- 



and as — 1°. 9. Using the revised, proper motion-based convergent point estimate, and the 
new estimate of the velocity dispersion of the group, has negligible effect on the distance 
determinations. For these reasons (and clarity of presentation), I have chosen not to list the 
reevaluated quantities. 

After excluding TWA 12, 17, 18, 19, 22, and 24 as TWA members, I characterize the 
final TWA membership with probability plots (so as to be immune to the effects of outliers; 
Lutz & Upgren 1980): d TWA = 49 ± 3 (s.e.m.) ± 12 (la) pc, a TWA = 174°. 8 ± 12°.3 (la), and 
&twa — -37°. 1 ± 7°. 5 (la). The projected radii in a and 8 correspond to ~7pc at d — 49pc. 
Taking into account the typical distance errors (~5 pc) and the observed distance dispersion 
(~12pc), the data are consistent with the radius along the line of sight being ~10pc (~40% 
larger than the projected width of ~7pc). All three of the new brown dwarf members (TWA 
26-28; red open circles in Fig. 3) cluster lie between d = 40-53 pc, close to the classic TWA 
membership. 

With the membership and characteristics of the TWA better defined, one can ask the 
question: are there other stars in the vicinity whose astrometric data also suggest that they 
are TWA members? Dozens of other young, low-mass, field stars have been proposed as 
TWA members, enough so that assessing their membership is probably worth a separate 
study. A question that can be answered here is: are there any high-mass TWA members 
besides HR 4796?. The quoted magnitude limits of the Hipparcos catalog suggest that it 
should be complete for unreddened A and B-type stars on, or above, the main sequence 
within ~85pc (Perryman et al. 1997). 1 queried the Hipparcos database for stars within 
a 15° radius centered on the TWA central position given earlier. I retained the 31 stars 
with parallaxes of >10mas and B — V colors of <0.30 (consistent with unreddened stars 
earlier than F0). I calculated membership probabilities and predicted cluster parallaxes for 
these stars, in the same manner that was done for 2M1207 in §2.4. Of these 31 stars, only 
eight had membership probabilities of >5%. For these eight, I compared the moving cluster 
parallax values to the Hipparcos trigonometric parallaxes. Only three of these eight stars 
had agreement between cluster and trigonometric parallaxes at better than 2a: HR 4796 
(known member), HIP 54477 (A1V, d = 58pc), and HIP 53484 (FOV, rf = 97pc). HIP 53484 
is ~4<t more distant than the mean TWA distance, and nearly ~15° from the TWA centroid 
position, so I reject its TWA membership, and discuss it no further. HIP 54477 is not so 
easy to dismiss as a TWA member. This Al dwarf has a high TWA membership probability 
(90%), and its trigonometric parallax (17.2 ± 0.7mas) agrees fairly well with its predicted 
TWA cluster parallax (20.8 ±1.6 mas). Its projected position is in the core region near TWA 
2, 4, and 8. At d ~ 56 pc (Hipparcos) or d ~ 48 pc (predicted cluster parallax distance), it 
would be slightly further than TWA 2, 4, and 8 (all of which have d ~ 40 pc). The radial 
velocity of HIP 54477 is not well constrained (v ra d = ±16.2 ± 10 kms -1 ; Barbier-Brossat & 



- 11 - 



Figon 2000), but consistent with that for a TWA member at its position (+12.2 ± 1.6kms _1 ). 
The star appears to be close to the zero-age main sequence, and so could be as young as 
the other TWA members. Further observations should be undertaken to see if the object 
has any low-mass companions which may further constrain its age. The membership of HIP 
54477 to TWA can not be rejected on kinematic grounds, but one would certainly like to see 
further data before claiming that it is a true TWA member. In summary, the TWA appears 
to contain at least one (HR 4796), but possibly two (HIP 54477), stars hotter than F0 in its 
membership. 



2.6. Expansion Age of TWA 

One may be able to put an interesting astrophysical constraint on the age of the 2M1207 
system through calculating an "expansion age" for the TW Hya association. MF01 claimed 
that the kinematics of the TWA are consistent with an expansion age of 8.3 Myr. The 
analysis of MF01 included tens of X-ray-selected stars in their analysis which have been 
since shown to not be pre-MS stars (SZB03, Torres et al. 2003). As the majority of the stars 
in the MF01 analysis are not genetically related to TW Hya or its cohort, this expansion age 
is not a useful constraint on the age of the TWA or 2M1207. With the best proper motion 
and radial velocity data currently available, I investigate whether an expansion is still evident 
in the TWA using a Blaauw expansion model. For discussions on trying to detect the linear 
expansions of unbound associations, see Blaauw (1956, 1964), Bertiau (1958), Jones (1971), 
Brown et al. (1997), Dravins et al. (1999), and Madsen et al. (2002). 



2.6.1. The Blaauw Linear Expansion Model 

Linear expansion of an association can not be demonstrated with proper motions alone 
(Blaauw 1964; Brown et al. 1997). A group of stars with generally parallel motion vectors, but 
with a small linear expansion, will simply appear to converge to a point further away (higher 
A) than that demonstrated by a group with strictly parallel motion vectors. The classical 
convergent point method equations which assume parallel motion are slightly modified to 
allow for expansion. In the Blaauw (1964) linear expansion model, the individual cluster 
parallax (w) for an association member is calculated as: 

\i v A 

and the radial velocity is predicted to follow the relation: 

Vrad = v' COS A' + K d + K (2) 
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where A is the AU as previously defined, fi v is the proper motion directed toward the 
convergent point, k is the expansion term in units of kms _1 pc _1 , d is the distance to the 
star in pc (where d pc ~ 1000 w^J, and K is a zero-point term which may reflect gravitational 
redshift or convective blueshift terms (see Madsen et al. 2003, for a detailed discussion). The 
"expansion age" r of the association in Myr is: 

T = 1 - l K- 1 (3) 

where 7 is the conversion factor 1.0227 pc Myr -1 km -1 s. The cluster speed v' and star-vertex 
angular separation A' are defined differently that in the standard case of parallel motions. 
In the Blaauw model, v' is the barycentric speed of a hypothetical association member 
participating in the expansion, situated at the barycenter of our solar system (see Fig. 3 of 
Blaauw 1964), and A' is the angular separation between a star and the association convergent 
point as defined solely by the stars' proper motions. If an association is expanding, the 
convergent point determined from the mean 3D space motion of its members (the "centroid" 
space motion) will define a different "convergent point" than the vertex determined through 
a convergent point analysis of the stars' proper motions. 

To test whether the association is expanding or not, and possibly assign an "expansion 
age", I analyze the data for TWA members two ways. First, I compare model convergent 
points for varying expansion ages to the observed convergent point. Second, I will use the 
available radial velocities and cluster parallax distances to directly measure the expansion 
rate. 



2.6.2. Expanding versus Non-expanding Association Convergent Point 

In Fig. 2, I plot the variation in the convergent point (long-dashed line) if one takes the 
TWA "centroid" space motion vector (using the mean velocity vector for TWA 1, 4, 9, and 
11), and add linear expansion with characteristic expansion timescales. In §2.5, I determined 
that the best convergent point for the final TWA membership using the proper motion data 
alone was (a = 100°. 5 ± 4°. 2, 5 = —27°. 9 ± 1°.9). Predicted expansion model convergent 
points for ages 0-100 Myr were statistically compared to the observed convergent point error 
ellipse. From this analysis alone, one can reject expansion ages of <7 Myr at 5% significance, 
and <6Myr at 1% significance. 

The close agreement between the TWA vertex found by the convergent point method 
and the vertices for the four individual TWA members with known UVW vectors (see Fig. 
2), suggests that any kinematic expansion must be very subtle, and perhaps not even demon- 
strable with existing data. It is worth exploring whether the radial velocity data can help 
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either determine a significant expansion age, or at least place a more interesting lower limit. 



2.6.3. A "Hubble Diagram" for TWA? 

Blaauw (1956, 1964) suggested that linear expansion or contraction may be detectable 
if deviations are present between the observed spectroscopic radial velocities, and those 
predicted from the moving group method for parallel motion. If a significant linear expansion 
term n is present, then the Blaauw expansion model equations (Eqns. 1 and 2) predict that 
one should see a correlation between distance d and the difference (v ra d - v cosA) between the 
observed and predicted spectroscopic radial velocities. As the radius of the TWA is ~10pc 
and the isochronal age is ~10 7 yr, one expects that k should be of order ~0.1 kms -1 pc" 1 , if 
the stars are linearly expanding from a point. 

The effects of expansion on cluster parallax distances are usually negligible (e.g. as 
shown by comparing trigonometric parallaxes to cluster parallaxes; e.g. de Bruijne 1999b; 
Mamajek et al. 2002; Madsen et al. 2002). For the case of TWA, the change in cluster 
parallax distances, between assuming parallel motion and linear expansion, is <6% rms for 
expansion ages of >5Myr, and <3% rms for >10 Myr. Note that expansion ages of <6Myr 
were effectively ruled out in §2.6.2, and the typical distance errors from other sources (e.g. 
proper motions) are ~11%. One can then conclude that the effects of association expansion 
(if any) on the distances and distance errors quoted in this study are negligible. 

In order to detect any possible expansion by fitting the Blaauw model to the observa- 
tions, I adopt the convergent point defined solely using the proper motion data. I estimate 
v' for the four TWA members (TWA 1, 4, 9, 11) with trigonometric parallaxes through the 
equation: 

v' = (4) 
zd sin A' 

The mean value for the four TWA members is v' = 20.4 ±2.2 km s -1 . Already, one notices 
that v (21.3 ± 1.3kms _1 ) is indistinguishable from v', which is consistent with no expansion. 

In order to see whether a non-zero k coefficient is detectable, I plot in Fig. 4 the data 
in the format d versus (v ra d — v' cosX') so that one can solve for the slope k and intercept K: 

Vrad — v' COS A' = K d + K (5) 

Plotted in this form, any expansion will manifest itself as a significantly positive slope. The 
individual distance estimates for the expansion model are calculated as: 

1000 v' sin X' 

d » c = Afi v (6) 
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As seen in Fig. 4, it is a success of the kinematic model that the (v ra d — v' cosX') values 
are crowded near zero at all. Recall that the predicted radial velocity component (V cosX') 
is totally independent of any measured radial velocity data, i.e. they are solely dependent 
on the the convergent point position (via A'), and the trigonometric parallax distances and 
proper motions for TWA 1, 4, 9, and 11 (via v'). This agreement further strengthens the 
interpretation that the TWA constitutes a bona fide kinematic group. 

The errors in distance and velocity difference have some peculiarities worth mentioning. 
The distance errors tend to scale with distance, i.e. cr^ oc d. Secondly, the distances will all 
be affected systematically if the convergent point is in error. Finally, the linear fit of the 
data to equation #5 using weighting in both variables (using fitexy from Numerical Recipes; 
Press et al. 1992) gives an uncomfortably good fit (x 2 /^ = 7.8/19), presumably due to 
overestimated errors in either the observed radial velocities, group speed, or convergent 
point. This weighted fit finds k = +0.036 ±0.039 km s' 1 pc^ 1 and K = +1.07 ± 0.51 kms" 1 , 
but again, due to the very low x 2 , it is unclear how much to believe the errors. 

To avoid overinterpreting a derived slope k whose error bars may not be believable, 1 fit 
an unweighted, ordinary least-squares line to the data, with the distance d as the independent 
variable, and the velocity difference {v ra( i — v' cosX') as the dependent variable. I do this for 
the 19 TWA "final" members whose radial velocity errors are <2.5kms _1 . Since the sample 
is small, I use bootstrap and jackknife resampling to help determine the error in the derived 
slope (Feigelson & Babu 1992), although the agreement with the errors derived from the 
asymptotic formulae is excellent. The least-squares fit finds k = +0.049 ± 0.027kms~ 1 pc _1 
and K = +1.20 ± 0.36 kms -1 (evaluated at the mean distance). Although the sign of the 
slope is consistent with expansion, the correlation is very weak (Pearson r = 0.42 + 0.19). 
The basic result is unchanged whether all of the TWA members are retained, independent 
of radial velocity error, or if only the 8 TWA stars with radial velocity errors of <2kms~ 1 
are retained 4 . 

Although the slope k is small, one can state that it is positive at 95% confidence, i.e. 
that the data are consistent with some expansion. Unfortunately, the derived expansion 



4 A non-zero velocity offset K should not cause too much alarm. Part of the offset may be due to 
gravitational redshift, which for the typical ~10Myr-old TWA member with mass ~O.5M should be of 
order ^0.4kms _1 (Greenstein & Trimble 1967, using radii from the D'Antona & Mazzitclli 1997 tracks), 
compared to ~0.6kms _1 for the Sun. An unexplained radial velocity offset of 0.4 kms -1 appears to be 
present among low-mass Hyades members (Gunn et al. 1988), even after accounting for gravitational redshift. 
The offsets between measured "spectroscopic" radial velocities and "astrometric" radial velocities arc difficult 
to quantify, but should be more easily measurable for larger samples of stars with future astrometric missions 
(Dravins et al. 1999; Madsen et al. 2003). 
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age has very large errors, and is of limited utility: r = 7 -1 k~ x ~ 20+y 5 Myr. The prob- 
ability distribution function of k excludes expansion ages of <8.7Myr at 99% confidence, 
and <10.4Myr at 95% confidence. The confidence intervals on the expansion age are very 
wide: 13-43 Myr (68% CL) and 9.5-262 Myr (90% CL), with ~4% of the probability dis- 
tribution corresponding to contraction. The expansion age advocated by MF01 (8.3 Myr) 
can, however, be statistically rejected. It does not seem appropriate at this time to quote 
an unambiguous "expansion age" for the TW Hya association, but to quote the lower limit 
(> 104 Myr). 



3. DISCUSSION 

With an improved distance estimate, one can revise the absolute magnitude, luminosity, 
and inferred mass estimates for 2M1207A and B. The properties of 2M1207 A and B from 
the literature, and derived here, are listed in Table 3. Using the photometry from Chauvin 
et al. (2004) and revised distance estimate from the moving cluster method, the absolute 
magnitudes of 2M1207A and B are M^(A) = 8.32 ± 0.27 and M^(B) = 13.30 ± 0.29 mag, 
respectively. These are 0.6 mag fainter than one would derive using d — 70 pc (i.e. a factor 
of two intrinsically dimmer). I calculate luminosities using these absolute magnitudes, and 
the bolometric correction estimates of Golimowski et al. (2004). Using the constraints on 
luminosity and age (Chauvin et al. 2005), I interpolate masses from the non-gray evolutionary 
tracks of Burrows et al. (1997), the DUSTY tracks of Chabrier et al. (2000), and the COND 
tracks of Baraffe et al. (2003). For all three sets of evolutionary tracks, the masses of A and 
B cluster near ~21M Jup and ~34M Jup . Table 3 also lists the mass extrema from the la 
extrema in both luminosity and age (i.e. the low mass end is for the -la luminosity and 
age, and the high mass end is the +la luminosity and age). With the previous distance 
estimates (~70pc), Chauvin et al. (2004) estimates mass of 25Mj up and 5Mj up for A and 
B, respectively. For all three models, the inferred upper mass limit of 2M1207 B (~5-7Mj up ) 
is less than half of the deuterium-burning mass limit (~13Mj up ; Burrows et al. 1997), and 
less than half of the maximum mass of Doppler velocity planets (~15Mj Mp ; Marcy et al. 
2005). Hence, 2M1207 B could be considered a "planet" on the merits of its inferred mass. 

The angular separation of AB (778 mas) measured by Chauvin et al. (2004) translates 
into a projected physical separation of 41 ± 5 AU at the revised distance (similar to the semi- 
major axis of Pluto). If the observed separation is assumed to be equal to the semi-major 
axis, and one adopts masses of 3.5 and 21Mj np for A and B, then one naively predicts an 
orbital period of ~1700yr. The pair has a high mass ratio (q ~ 0.2), and B is massive 
enough to force the primary to be ~6 AU from the system barycenter. A solid detection of 
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orbital evolution (or any hint of the dynamical masses of the components) will probably not 
be reported anytime soon. 

One can not rule out whether the TWA is expanding on a timescale longer than its 
isochronal age (>10Myr). The slow, or negligible, expansion may also be a clue that the 
proto-TWA molecular cloud complex was perhaps not a small pc-sized core with tens of stars, 
similar to those seen in Taurus (e.g. LDN 1551). The TWA members may have formed in 
a series of small-N systems (N ~ few stars) distributed along filaments, separated by a few 
pc, and with similar bulk motions. The TWA appears to be moving away from the LCC 
subgroup (Mamajek et al. 2000), so it is conceivable that the proto-TWA cloudlets were 
simply fragments of the proto-LCC cloud, which owed their velocities to molecular cloud 
turbulence (Feigelson 1996). An alternative scenario is that the proto-TWA cloudlets were 
bright-rim clouds or cometary globules on the periphery of LCC ~10-15 Myr ago, when 
presumably the LCC subgroup still had a few late-0 stars (de Geus 1992). Such cloudlets 
could have been accelerated away from the LCC O-star population through the rocket effect 
(Oort & Spitzer 1955), and compressed to form stars due to radiation-driven implosion 
(Bertoldi & McKee 1990). The energy input from deceased LCC members (via UV light, 
winds, and supernovae) has probably dominated the energy input of the local interstellar 
medium over the past 10 Myr, and within 100 pc, in the general direction of LCC and TWA 
(Maiz-Apellaniz 2001). Small-scale star-formation in cometary globules on the edge of OB 
associations has strong observational support (Reipurth 1983; Ogura & Sugitani 1998), and 
there is strong evidence for triggering by the massive stellar population (e.g. Kim et al. 
2005; Lee et al. 2005). A cometary globule formation scenario for TWA might explain a few 
observational quirks of the group, namely its location (~70 pc away from the nearby LCC OB 
subgroup), age (~7 Myr younger than LCC), space motion vector (directed ~5kms~ 1 away 
from the LCC; Mamajek et al. 2000), and low stellar density. The small, young stellar groups 
associated with rj Cha, e Cha, and (3 Pic show many of these same symptoms (Mamajek et 
al. 1999, 2000; Ortega et al. 2002; Jilinski et al. 2005), although the r] and e Cha clusters 
appear to be more strongly bound than the TWA and (3 Pic groups. Cloudlets analogous to 
those on the periphery of Vel OB2 (Kim et al. 2005) and Ori OBI (Lee et al. 2005) may be 
the evolutionary predecessors of small, unbound, ~10 Myr-old associations like TWA. That 
2M1207 and the TWA formed in a region of rather low stellar density could explain how 
such a wide, low-mass binary system as 2M1207 could survive its birth environment intact. 
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Fig. 1. — The relative change in position of 2M1207 since the late 1970s in a (top) and 5 
(bottom). The position offsets are relative to the positions on the ICRS at mean epoch. The 
mean epoch positions and mean epochs are a Q = 181°. 889563 (t (a ) = 1995.08) and 5 Q = 
-39°. 548315 (t (5 ) = 1996.57). The effects of parallax were not included since the predicted 
amplitude (~20 mas) is smaller than the individual position errors. The positional errors at 
mean epoch are a ao ~ o~s ~ 50 mas. The data are consistent with /i Q * = —71.6 ± 6.7masyr~ 1 
and yu 5 = -22.1 ± 8.5 mas yr" 1 for 2M1207. 
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Fig. 2. — New and previously determined convergent point (vertex) estimates for the TW 
Hya association. The open circles are the inferred convergent points for TWA 1, 4, 9, 11, 
and 19 based on their UVW space motions (and la error bars). Large Xs are previously 
published TWA vertices from Frink (2001, F01), MF01, R03, and SZB03. Filled triangles 
are the vertices for several subsamples of TWA objects (TWA 1-13, 14-19, 21-25, and 14-25) 
found through the convergent point method (based solely on the proper motion data). The 
dotted line and dashed line are the 68.3% (la) and 95.5% (2a) confidence levels in a and 
6 around the vertex found for the "classic" TWA sample (TWA 1-13). The error regions 
around the other filled triangles (TWA 14-25, 21-25, 14-19) are 50-100% larger than that 
for TWA 1-13, and similarly shaped (but not shown for clarity). Predicted vertices for the 
TWA with a wide range of expansion ages are distributed along the long dashed line. 
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Fig. 3. — Right Ascension versus distance for candidate members of TWA and LCC. Solid 
black circles are TWA 1-11 from Webb et al. (1999b), solid black squares are TWA 12-13 from 
Sterzik et al. (1999), open green triangles are TWA 14-19 from Zuckerman et al. (2001), open 
cyan square is TWA 20 from Reid (2003), open blue pentagons are TWA 21-25 from Song et 
al. (2003), open red circles are TWA 26-28 (2M1139, 2M1207, SSSPM J1102), small Xs are 
pre-MS LCC members from Mamajek et al. (2002), and small open squares are B-type LCC 
members from de Zeeuw et al. (1999). The dotted line is the mean TWA distance relation 
from Reid (2003). The dashed line at d = 85 pc is shown to illustrate the detached nature of 
the TWA and LCC groups, i.e. although they have similar space motions, they do appear 
to occupy separate regions. If the likely LCC members (TWA 12, 17, 18, 19, and 24) are 
excluded, as well as the probable nonmember TWA 22, the rest of the TWA membership is 
consistent with having distances of <i = 49± 12 (la) pc. 




Fig. 4. — Distance versus the difference between the observed radial velocity (*w) and 
the distance-dependent radial velocity component of the Blaauw expansion model (v cosA) . 
"Final" TWA members are filled circles, and probable LCC members are open circles. For 
parallel motion and no expansion, a slope of zero is expected. In the case of expansion, one 
expects the closest members to be more blueshifted, and the more distant members to be 
more redshifted. The solid line is the best fit slope k (= 0.049 ± 0.027kms _1 ) to the data, 
with la error bars as dotted- dashed lines. The prediction for a system with expansion age of 
8.3 Myr (MF01) is shown by the dashed line, and ruled out by the data. 
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Table 1. Properties of Proposed TW Hya Association Members 







W 


1 A\ 
V 4 ) 


(5J 


(6 J 


V) 


(°) 


l y ) 


(,10) 


fi 1 *i 
( LL ) 


TWA 




fJ>a* 




rxCI. 




,i 
<i 


v rad(W e d-) 


v rad(°^> s -) 


rxei. 


r mai 


44- 


Name 


[masyr -1 ] 


[masyr -1 ] 




[/oj 


lP C J 


ri,™ — li 
[Km s j 


[Km s j 




Member ? 


1 


TW Hya 


-66.8 ±1.6 


-15.2 ±1.3 


1 


97 


51 ±4 


±13.2 


+12.7 + 0.2 


7 


Y 


2 


CD-29 8887 


-91.6 ±1.8 


-20.1 ± 1.3 


2 


57 


39 ±3 


+12.4 


+11.0 + 0.1 


7 


Y 


3 


Hen 3-600 


-109.3 ±8.7 


±0.8 ±8.9 


3 


18 


34 ±4 


±12.9 


±12.5 ±1: 


8,9 


Y 


4 


HD 98800 


-91.7 ±1.5 


-30.0 ±1.5 


1 


73 


40 ±3 


±11.0 


±9.3 ±1: 


10,11 


Y 


5AB 


CD-33 7795 


-85.4 ±3.6 


-23.3 ±3.7 


1 


88 


44 ±4 


+11.4 


var. 


7 


Y 


6 


TYC 7183-1477-1 


-57.0 ±2.1 


-20.9 ±2.1 


1 


53 


51 ±5 


+15.4 


±16.9±5 


9 


Y 


7 


TYC 7190-2111-1 


-122.2 ±2.3 


-29.3 ±2.3 


4 


100 


27 ±2 


+14.2 


+11.7±2 


9 


Y 


8A 


GSC 06659-01080 


-99.3 ±9.0 


-31.3 ±8.9 


3 


83 


38 ±4 


+10.5 


±7.8 ±2 


9 


Y 


8B 


2MASS J11324116-2652090 


-95.3 ±10.0 


-29.5 ±10.3 


3 


84 


39 ±5 


+10.5 


±7.8 ±2 


9 


Y 


9A 


CD-36 7429A 


-52.8 ±1.3 


-20.2 ±1.8 


1 


81 


69 ±6 


+10.6 


±9.5±0.4 


7 


Y 


9B 


CD-36 7429B 


-70.7 ±13.3 


-6.6 ± 15.8 


3 


67 


56 ±12 


+10.6 


+11.3 ±2 


9 


Y 


10 


GSC 07766-00743 


-72.6 ± 12.2 


-32.1 ±12.3 


3 


100 


53 ±9 


+8.3 


+6.6 + 2 


9 


Y 


11A 


HR 4976A 


-53.3 ±1.3 


-21.2 ±1.1 


4 


99 


73 ±6 


+8.0 


+6.9 + 1.0 


9,12,13 


Y 


12 


RX J1121. 1-3845 


-36.3 ±8.6 


-1.6 ±8.9 


3 


69 


103 ±26: 


+12.4 


+10.9 + 1.0 


7,14 


N 


13 


RX J1121. 3-3447 


-67.4 ± 11.8 


-17.0 ± 11.8 


3 


99 


53 ±10 


+12.0 


+12.1 ± 1: 


7,14 


Y 


14 


UCAC2 12427553 


-43.4 ±2.6 


-7.0 ±2.4 


1 


96 


80 ±8 


+13.1 


+16.0 + 2 


9 


Y? 


15 


GSC 08236-01074 


-100.0 ±33.0 


-16.0 ±6.0 


5 


88 


41 ±6 


+9.1 


+11.2 + 2 


9 


Y 


16 


UCAC2 12217020 


-53.3 ±5.2 


-19.0 ±5.2 


1 


100 


72 ±9 


+8.8 


±9.0 ±2 


9 


Y 


17 


GSC 08248-00700 


-28.0 ±8.5 


-11.1±8.5 


3 


84 


163 ±46: 


+6.3 


±4.6 ±6 


9 


N 


18 


UCAC2 12908626 


-29.0 ±5.2 


-21.2 ±5.2 


1 


64 


121 ±20: 


+6.0 


±6.9 ±3 


9 


N 


19A 


HD 102458A 


-33.6 ±0.9 


-8.5 ±0.9 


1 


73 


109 ±9: 


+11.6 


+13.5 + 2.4 


7,9 


N 


19B 


HD 102458B 


-35.6 ±4.8 


-7.5 ±4.6 


1 


97 


103 ±16: 


+11.6 


+15.2 + 2 


9 


N 


20 


GSC 08231-02642 


-52.0 ±5.0 


-16.0 ±6.0 


5 


97 


75 ±9 


+9.0 


±8.1 ±4 


15 


Y 


21 


HD 298936 


-65.3 ±2.4 


±13.7±1.0 


1 


99 


45 ±4 


+15.8 


+17.5 + 0.8 


16 


Y 


22 


SSSPM J1017-5354 


-176.0 ±7.0 


-22.0 ±8.0 


6 


2 


18 ±2: 


+15.5 






N? 


23 


SSSPM J1207-3247 


-68.0 ±4.0 


-23.0±4.0 


6 


86 


57 ±5 


+8.9 






Y 


24 


MML 5 


-34.4 ±2.8 


-13.1 ±1.7 


1 


10 


107 ±12: 


+11.1 


+11.9 + 0.9 


16 


N 


25 


TYC 7760-283-1 


-75.0 ±2.0 


-26.9 ±1.4 


1 


100 


51 ±4 


+9.2 


+9.2 + 2.1 


16 


Y 


26 


2MASSW J1139511-315921 


-93.0 ±5.0 


-31.0 ±10.0 


6 


99 


40 ±4 


+10.6 


+11.6 + 2 


17 


Y 


27 


2MASSW J1207334-393254 


-71.6 ±6.7 


-22.1 ±8.5 


3 


98 


53 ±6 


+9.7 


+11.2 + 2 


17 


Y 


28 


SSSPM J1102-3431 


-82.0 ±12.0 


-12.0 ±6.0 


6 


70 


43 ±7 


+13.2 






Y 



Note. — Columns: (1) TWA number, (2) other name, (3) proper motion in RA (/i a » = fi a cos<5; ICRS frame; 1<t errors), (4) 
proper motion in Dec (fi^; ICRS frame; la errors), (5) proper motion reference, (6) membership probability (§2.4), (7) predicted 
distance from moving group method (§2.5) and la uncertainty, (8) predicted radial velocity from moving group method (with 
uniform 1.6 km s^lcr uncertainty), (9) observed mean radial velocity, (10) radial velocity reference, (11) final TWA membership 
assessment. For some binaries, I have estimated the systemic velocity by assuming a mass ratio. These are probably good to 
~1 kms , and their errors have been marked with colons. For the systemic RV of HD 98800, I adopted the component masses 
from Prato et al. (2001) and assumed M(Ab) = 0.5 Mq. The radial velocity listed for TWA 15AB is that measured for A, but 
the proper motion is for the photoccntcr of AB (TWA 15A = 2MASS J12342064-4815135, TWA 15B = 2MASS J12342047- 
4815195). Kinematic distances to non-members (TWA 12, 17, 18, 19AB, and 24 are probably LCC members; TWA 22 may 
not be a member) should not be taken seriously. However, due to the similarity in space motions between TWA and LCC, the 
distances to TWA 12,17,18,19AB, and 24 listed are probably within 5% of the real distance if they are indeed LCC members. 
References: (1) Zacharias et al. (2004, UCAC2), (2) Platais et al. (1998, SPM), (3) calculated by the author using positions 
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from the following catalogs: GSC-ACT (Lasker et al. 1999), USNO-A2.0 (Monet et al. 1998), GSC 2.2 (STSci & Osservatorio 
Astronomico di Torino 2001), 2MASS (Cutri ct al. 2003), DENIS (The DENIS Consortium 2003), UCAC1 (Zacharias ct al. 
2001), and Rousseau et al. (1996), (4) Tycho-2 (H0g et al. 2000), (5) USNO-B1.0 (Monet et al. 2003), (6) Scholz ct al. (2005), 
(7) Torres et al. (2003), (8) dc La Reza et al. (1989), (9) Reid (2003), (10) Torres et al. (1995), (11) Prato ct al. (2001), (12) 
Barbier-Brossat & Figon (2000), (13) Grcnier et al. (1999), (14) Stcrzik ct al. (1999), (15) Webb (1999, as reported in Reid 
(2003)), (16) Song ct al. (2003), (17) Mohanty ct al. (2003) 
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Table 2. Trigonometric vs. Cluster Parallaxes 



(1) (2) (3) 

Name ^trig ^clus 

[masl fmasl 



TW Hya 
HD 98800 
TWA 9 
HR 4796 



17.8±2.2 
20.5±2.8 
19.9±2.4 
15.1±0.7 



19.5 ±1.6 
25.0 ±2.0 
14.8±l.l a 
13.8±1.1 



Note. — Columns: (1) com- 
mon star name, (2) weighted 
mean of Hipparcos and Tycho- 
1 trigonometric parallax (Perry- 
man et al. 1997), (3) parallax 
from cluster parallax method 
(this work). 



a Weightcd mean of individual 
estimates for TWA 9 A and B. 
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Table 3. Properties of 2M1207 A and B 



(1) 


(2) 


(3) 


(4) 


(5) 


Property 


2M1207 A 


2M1207 B 


Units 


Notes 


Spectral Type 


M8.5±l 


L7.25±2.25 




1 


K 


11.96 ±0.03 


16.93±0.11 


mag 


2 


M K 


8.32 ±0.27 


13.30 ±0.29 


mag 


3 


BC K 


3.12±0.14 


3.25 ±0.14 


mag 


4 


log(L/L ) 


-2.68 ±0.12 


-4.72 ±0.14 


dcx 


5 


Mass (Baraffe et al. 2003) 


21(19-30) 


3.3(2.3-4.8) 


M Jup 


6 


Mass (Chabrier et al. 2000) 


21(19-31) 


3.2(2.3-4.8) 


Mj„ p 


7 


Mass (Burrows et al. 1997) 


20(17-24) 


4.2(2.6-6.5) 


M Jup 


8 



Note. — (1) spectral types from Chauvin et al. (2004), (2) iiT-band 
photometry from Chauvin et al. (2004), (3) absolute K magnitudes 
assuming distance from §2.5, and no reddening, (4) bolometric correc- 
tions are from polynomial of Golimowski et al. (2004). Error includes 
uncertainty in spectral type and rms of their BCx(SpT) fit, (5) lu- 
minosity using Mks and BCk and assuming Mj, ;Q = 4.75, (6) mass 
(and mass range) from COND evolutionary tracks of Baraffe et al. 
(2003). (For all evolutionary tracks, the best interpolated mass es- 
timate at age 8 Myr is given, followed by the extrema of the mass 
range considering the uncertainties in luminosity and age, where I 
assume an age of 8^3 Myr, following Chauvin et al. 2004), (7) mass 
(and mass range) from DUSTY evolutionary tracks of Chabrier et al. 
(2000), (8) mass from evolutionary tracks of Burrows et al. (1997). 



